Unified description of bulk and interface-enhanced spin pumping 
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The dynamics of non-equilibrium spin accumulation generated in metals or semiconductors by rf 
magnetic field pumping is treated within a diffusive picture. The dc spin accumulation produced 
in a uniform system by a rotating applied magnetic field or by a precessing magnetization of a weak 
ferromagnet is in general given by a (small) fraction of hui, where to is the rotation or precession 
frequency. With the addition of a neighboring, field-free region and allowing for the diffusion of 
spins, the spin accumulation is dramatically enhanced at the interface, saturating at the universal 
value hu) in the limit of long spi n relaxa tion time. This effect can be maximized when the system 
dimensions are of the order of ^2nD/u>, where D is the diffusion constant. We compare our results 
to the interface spin pumping theory of A. Brataas et al. [Phys. Rev. B 66, 060404(R) (2002)]. 
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A new direction in the field of spin electronics is 
the study of the interaction between spin currents and 
the magnetization dynamics of a ferromagnetic electrode 
into, or out of which the spin currents flow. In spin 
pumping devices 0,00; electron spins are pumped from 
a ferromagnet in resonance with an rf magnetic field, into 
an adjoining non-magnetic metal without accompanying 
charge currents. The so-called spin battery is based 
on the interface spin-mixing conductance which de- 
scribes the coherent reflection and spin rotation within an 
exchange length (typically a few nm in conventional fer- 
romagnets) of the interface. It has been experimentally 
verified by measuring the effect of adjoining metal lay- 
ers on the Gilbert damping constant of a ferromagnetic 
layer 0, [||, due to the transfer of angular momentum 
away from the ferromagnetic and into adjacent metals 
7J. For the spin battery, under specific conditions the 
universal value of fttu is obtained for the spin accumula- 
tion, where u> is the precession frequency of the magne- 
tization. On the other hand, a bulk-oriented approach 
based on mean-field theory [8( has been used to relate 
the intrinsic Gilbert damping to the spin relaxation pro- 
cesses resulting from the generation of non-equilibrium 
spin accumulation in the ferromagnet. 

We present a model in which spin pumping with a ro- 
tating magnetic field can produce spin accumulation in 
bulk material. Moreover, in a hybrid system in which 
spins are allowed to diffuse from a " pumped" region into a 
neighboring region without fields, the spin accumulation 
can be dramatically enhanced near the interface. In both 
cases the spin accumulation reaches the universal limit 
Hlo for long spin relaxation times. Our approach is based 
on classical dynamics of spin ensembles in a diffusive sys- 
tem, and is applicable to spin pumping in metals and 
semiconductors driven by a rotating, applied magnetic 
field, as well as in weak ferromagnets driven by magneti- 
zation precession. However, we emphasize that our semi- 
classical approach cannot be applied directly to the con- 



ventional strong ferromagnet regime in which the large 
exchange fields require a non-diffusive, quantum mechan- 
ical treatment because of the very small magnetic co- 
herence length, much less than the mean- free-path. The 
dc spin-accumulation induced by spin-pumping is first 
derived for a bulk, uniform system, for which we ob- 
tain the surprising result that bulk spin accumulation 
is generated even for small Zeeman splitting. Then, by 
connecting a magnetic field-free region and solving the 
boundary-value problem of diffusion across the interface, 
we find a transition from bulk to interface-enhanced spin 
accumulation. 

We start with Bloch-type equations written in terms of 
the chemical potentials for the three spin directions, / = 
(f x , f y , f z ) that describe the space- and time-dependent, 
non-equilibrium spin accumulation f(x, t) (restricting 
the description to one spatial dimension x) The 
basis of our description is the following equation for the 
dynamics of f(x,t) in a rotating magnetic field B(x,t) 

mm. 



dt 



T{x, t) = -DV 2 f+ f - + x f, (1) 



where D is the diffusion constant and I(x,t) is a source 
term which will be described in the next paragraph. The 
length scales relevant to the problem are the spin dif- 
fusion length A = \J Dt, the Larmor precession length 
Ab = \j2mD I\jJb\i where HCSb = 9HbB, and a length 
\ u = y/2wD/u) corresponding to the distance the spins 
can diffuse in one period — . Eqn. ^ with the left 
hand side equal to zero, has been shown to give an accu- 
rate description of spin dynamics in metallic systems for 
time-independent magnetic fields 0, 0] . 

The key new ingredient of our approach is the source 
term: 



I{x,t) = ~gfi B g^B(x,t). 



(2) 
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FIG. 1: (a) Magnetic field vector diagram. An applied rf 
field B xy rotates with angular frequency u in the x-y plane. 
A static magnetic field B z is applied perpendicular to B xy 
generating a field vector B that rotates around B z with cone 
angle 6. (b) The general hybrid system with an interface at 
x = 0. The fields just described are applied in the "pumped" 
region I, x < 0. B — in region II, x > 0. 



It describes the rate at which the oscillating magnetic 
field, via the Zeeman energy, pushes spins aligned with 
the magnetic field below the Fermi level, and anti-aligned 
spins above from the Fermi level. This mechanism can 
be seen as a source of locally-injected, time-dependent 
spin currents 0|. Without this source term, / = in 
the steady state. 

The role of Eqn. [3 is sketched out as follows: consider 
two magnetic fields applied in orthogonal directions, say 
B x and B y , equal in magnitude but oscillating out of 
phase by 7r/2: B x — B xy cosujt and B y = B xy smujt. 
The resulting I x generates spin accumulation in f x in 
phase with B y , and subsequent precession yields spin ac- 
cumulation in f z . Likewise, I y produces an f y that is ir 
out of phase with B x , and subsequent precession yields 
spin accumulation in f z with the same sign as for the 
I x case. Both f x and f y oscillate about a mean of zero 
with frequency uj. However, f z acquires a finite dc com- 
ponent. The steady-state solution of / can be imagined 
as a spin accumulation vector that follows the rotating 
field vector, lagging by phase angles which depend on ui 
and r. 

Fig. \Ql shows the basic field configuration we will 
consider: a rotating magnetic field B xy and a per- 
pendicular, static field B z . The field vector is B = 
(B xy cosut, B xy sinwt, B z ) in Cartesian coordinates. We 
first consider region I (Fig. \\]p) alone; i.e., a uniform 
(spatially invariant) system without boundaries. Drop- 
ping the diffusion term and rewriting Eqn^ in a more 
concise form: 



t/ = -/ - uj b t X / - Hu) b t, 



(3) 



where uj b = (uj xy cosujt,uj xy smujt,uj z ). Since the in- 
duced spin accumulations f x and f y will, in the steady 
state, oscillate with the same frequency as B x and B y , 
it will be advantageous to transform to a coordinate sys- 
tem in which the x and y axes rotate about the z axis 
with angular frequency uj = ujz, so that B = (B xy , 0, B z ). 



The time derivative of a vector r in the stationary frame 
transforms as (r)i a b =^ (r)rot + & x {r)rot- Using this 
relation, EqnOlnow becomes 



rf 



-f- (uj B + &)t x / - h(uj B +LU x uj b )t, (4) 



In the last term, we have uj b = in this frame and Q x 
lu b = (0, uitOxy, 0). For the stationary analysis we set 

/ = 0, and obtain the result 0: 



1 + {uJxyT) 2 + ((UJ Z + Lu)t) [ 



rfajj. 



(5) 



The other components are related to f z : f x = — ^^ + ") f z 
and f y = —(uj xy T)~ 1 f z . Note that, although we will not 
discuss these components further in this letter, they are 
an integral part of the model and crucial for the deter- 
mination of f z . An expression similar to Eqn. [S] has 
been used to describe Gilbert damping in a ferromag- 
net @, and Eqn. can be viewed as a generalization 
of that result. In resonance when uj z = —uj, the uni- 
versal value hw is achieved when uj xv t » 1. Here- 
after, we will focus on the more general case, where a 
finite (w + lu z ) strongly damps the spin accumulation 
|l4|. Then, in the limit of long spin relaxation times 
such that t -1 << (ui xy , ui z + uj), the prefactor reduces to 
(1 + (^^fr^) 2 ) -1 and is independent of r; f z will then at- 
tain the universal value of hw only when uj xy > > (uj z +uj) . 

We now turn to discussing the hybrid system of Fig. 
The rotating field is only applied in the space 
x < (the "pumped" region I), whereas for x > 
there are no magnetic fields (region II) . We will consider 
two physical realizations for the pumped region I: spin 
pumping in a non-magnetic material with an applied, ro- 
tating magnetic field of magnitude B xy = 10 mT, and 
B z = (both region I and II are non-magnetic conduc- 
tors); and a model of a weak ferromagnet with B xy = 1 T 
and B z = 100 T representing the exchange interaction 
corresponding to small angle magnetization precession 
around an easy axis (region II is a non-magnetic conduc- 
tor). This exchange interaction approach, in which the 
magnetization of the localized d-moments couples to the 
s-electrons and polarizes them, has been used in early 
spin-pumping theory |l5| (which was applied to trans- 
mission electron spin resonance), and more recently in 
the mean-field-theory description of Gilbert damping |]| . 

The system is solved as a boundary- value problem gov- 
erned by Eqn. 0]with the diffusion term reinstated: 



A 2 V 2 / = / + (uj b + uj)t x / + Huj x UJ B T. 



(6) 



This equation can either be solved numerically directly, 
for instance using a finite element method solver, or by 
using the boundary conditions to determine the six con- 
stants of integration (for each region) of the general so- 
lution 16]. For simplicity in our model system we have 
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FIG. 2: The spin accumulation f z /hu as a function of x/A 
away from the interface, for r ranging in decades from 100 ps 
to 10 s. Both regions are unbounded. For (a) , B xy — 10 mT 
and B z = 0, for (b) B xy = 1 T and B z = 100 T. 



assumed that the spin diffusion lengths, times and dif- 
fusion constants are the same for both regions. With 
these assumptions the boundary conditions at the inter- 
face reduce to continuity of / and V/ at x — 0. When 
the regions are unbounded, the outer boundary condi- 
tions are that in region I, the bulk result (Eqn. [3] and 
below) is recovered; in region II, f(x >> 0) = 0. When 
the regions are bounded, the boundary condition is that 
there is no leakage of spin current, so V/ = at the 
boundaries. In the calculations the diffusion constant 
has been fixed at D = 0.002 m 2 / s (such as for Al) so the 
spin diffusion length A scales as y/r, and uj = 10 GHz. 

We first consider the case where both regions are un- 
bounded. Fig. shows the results of the calculation 
for the two different field pumping regimes, plotting the 
spatial dependence of f z /fko near the interface for val- 
ues of r ranging over many decades. The most striking 
feature in Fig. [5] is that the spin accumulation at the 
interface is strongly enhanced over the bulk, uniform sys- 
tem value (at x « 0), and in the limit of very large r 
achieves the universal value ftu>, similar to the results of 
Ref. . The calculations seem to describe a spin current 
injected at the interface which diffuses into both regions, 
even though the pumping occurs uniformly throughout 
region I. The enhancement of the spin accumulation at 
the interface is the primary result of this letter. In the 
unbounded system this effect is universal only for un- 
physically large r, so we now try bounding the system 
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FIG. 3: Spin accumulation at the interface for the symmetri- 
cally bounded system plotted as a function of the normalized 
bounding length L/X^, where A„ = yj 2nD /u ~ 1.1 fim, for 
r ranging in decades from 10 -10 (bottom curve) to 10 s (top 
curve) . 



in order to optimize the effect. In what follows we will 
focus only on the weak ferromagnet model, however most 
of the results apply generally to both systems. Fig. |3| 
shows the spin accumulation at the interface obtained 
for symmetric bounding of the regions at x = ±L as a 
function of L/X u , over decades of r from 10~ 10 to 10 
s. The spin accumulation is optimized when the sample 
dimension L ~ 0.6A W , independent of r. Fig. 01 shows 
the dependence of the spin accumulation on the bound- 
ing length, calculated under various bounding conditions 
for a particular r = 10~ 7 s. Bounding only region II 
yields optimal spin accumulation for L/A w close to one, 
however upon bounding only the pumped region the spin 
accumulation peaks at a smaller length scale, correspond- 
ing to L ~ Xb- Bounding region II at L2 = 0.6A W and 
allowing L\ to vary yields a peak at L\ ~ 0.6Ab. The 
spin accumulation obtained for different specific bound- 
ing lengths as a function of r is summarized in Fig. 
The main message from this figure is that by engineering 
the sample dimensions it is possible to achieve universal- 
ity at much lower r. 

A qualitative explanation for the interface enhance- 
ment is the averaging effect region II has on the x- and 
y-components of the spin accumulation. These com- 
ponents diffuse into region II but can diffuse back after 
some delay time. The back-flow will have a distribution 
of spin angles with respect to the rotating field in region 
I that tends to average out their influence, an effect sim- 
ilar to increasing the spin relaxation rate t^ 1 for the x- 
and y-components of the spin accumulation. However 
f z does not experience such an averaging effect. In a 
simple model of spin accumulation in a uniform system 
it can be shown that reducing r xy while staying in the 
limit uj^,T xy T z >> 1 leads to an enhancement of the spin 
accumulation 161. 
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FIG. 4: Spin accumulation at the interface vs the normalized 
bounding length L/X w , calculated with r = 1CP 7 s under 
various bounding conditions indicated as L\,L2 in the figure. 
The largest enhancement is obtained for the L, O.GXuj case 
with L ~ 0.6A S = 0.6 V^ttD/ojb = 15 nm. 
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E, as shown for the fit to curve E (with 10 -5 as the fit- 
ting parameter for 77), even though our model is derived 
for an entirely different regime. 

Our model of spin pumping in a non-magnetic region 
by a rotating magnetic field can be applied quite gener- 
ally to diffusive systems, such as metals and semiconduc- 
tors. Considering a semiconductor with t as large as 
100 ns and experimentally quite accessible oj = 10 12 Hz, 
then a spin accumulation close to the universal result of 
fkv corresponding to 700 /zV can be achieved even in the 
non-resonant case. The measurement of small spin accu- 
mulation signals due to the spin Hall effect in semicon- 
ductors by Kerr rotation microscopy has recently been 
demonstrated |l7|. The spin pumping- induced spin ac- 
cumulation and its interface enhancement that we have 
described should be measurable in similar systems. 

We acknowledge useful discussions with G. E. W. 
Bauer and M. Costache. Support has been provided 
by the Stichting Funtamenteel Onderzoek der Materie 
(FOM) and the RTN Spintronics Network. 




FIG. 5: The spin accumulation at the interface vs r for spe- 
cific bounding lengths indicated as L\ , L2 (curves A through 
E). The open circles are based on Ref. with r\ = 

10" 5 ^ tanh (see text for details). 



In the interface scattering matrix model Q, the spin 
accumulation close to the interface is given as f z = 
hw j^rgjp- , where 9 is the precession cone angle and 
rj = (Tj/r) tanh(L/A)/(L/A) is a reduction factor which 
depends on the normal metal properties and the spin- 
injection rate Tj which is inversely proportional to the 
mixing conductance g^. In the limit of very large r and 
A, the reduction factor goes to zero and the cone angle 
dependence factors out, leaving the maximum spin accu- 
mulation f z = fku. The authors restrict their analysis 
to the case where the ferromagnet thickness should be 
less than the spin diffusion length. To compare to their 
result we take 9 ~ -?p- = yi^ and evaluate t\ for L = A w , 



tanh 



This 



and find that ij depends on r as r\ cx -4= 
functional dependence on r agrees well with curves D and 
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